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We re-examine Quantum Electrodynamics (QED) with massless electron as a finite quantum 
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1 Introduction 



Subsequent to the work of Gell-Mann and Low ||l| in their classic 1954 paper, the possibihty that 
spin-| Quantum Electrodynamics (QED) may be a finite quantum field theory has been investigated 
by Johnson, Baker and Wiley and by Adler 0|. The basic premise of such a finite field theory 
is that the renormalization constants of QED, to wit, the electron bare mass ttlq, the electron wave 
function renormalization constant Z2 and the photon wavefunction renormalization constant Z3 
tend to finite limits as the cut-off A used in the calculations of the theory tends to infinity. It has 
been shown by Johnson et al that the renormalization constant Z2, a gauge variant quantity, can 
be rendered finite by an appropriate choice of gauge, while the gauge invariant quantity mo obeys 
the scaling law 

/ A2 \ "' 3 3 

mo (A) = const. ^ , e = -(ao/27r) + -(ao/27r)2 + ■ ■ • (1) 

y m^ y 2 8 

where ao is the bare fine-structure constant in QED, realized as an eigenvalue of the Gell-Mann- 
Low equation, il){x) = at Xq = ao- The bare mass ttlq thus vanishes as A — >■ 00, provided e > 0. 
This implies that the physical mass m arises solely from its interaction. With m-o = 0, as A — > 00, 
the axial-vector current J^ = ip'jfx'y^ip is formally conserved, except for the Bell-Jackiw-Adler |Q 
anomaly: d'^J^ = 2imoJ^ = 0. It was shown by Baker and Johnson [^ in 1971 that although 
the bare mass tends to zero in the limit of infinite cut-off, the matrix elements of J^ diverge in 
just such a fashion as to render the matrix elements of (9^ J^ nonvanishing and finite. In this case 
therefore the chiral symmetry is explicitly broken in spite of the vanishing bare mass. Hence there 
is no Nambu-Goldstone boson and such a zero mass state is not observed in nature either. 

To demonstrate that one can construct a finite theory of QED that will allow non-trivial eigen- 
values of the Gell-Mann-Low equation has remained somewhat of an outstanding problem. In our 
earlier work on the subject 0, we had concluded that such a self-consistent finite theory does not 
exist, reinforcing the 1979 conjecture of Baker and Johnson 0, ^. On the other hand there have 
been other investigations which conjecture that such a theory is feasible, especially the existence 



of a strong coupling phase of QED [^ |T^, |TT|]. We thought it worthwhile to re-examine the ques- 
tion of non-trivial solutions of the Gell-Mann-Low function in the finite QED theory, and embark 
on another careful investigation of the standard finite QED on the basis of the Dyson- Schwinger 
equation. 

In this investigation we shall begin with a solution to the Ward-Takahashi (WT) identity satisfied 
by the vector vertex function in QED and develop the Dyson-Schwinger equation satisfied by the 



electron propagator by employing the gauge technique. This will be presented in Sec. 2. We shall 
circumvent the problems introduced by the Bell- Jackiw-Adler anomaly by not dealing with the axial- 
vector current. We shall be working in the Landau gauge throughout. The solution to the WT 
identity for the vector vertex function in QED is arbitrary upto the existence of an undetermined 
transverse part, as is well-known. We shall not neglect such a transverse part and we shall see 
that it plays a crucial role in our investigation. In Sec. 3, we shall present a complete and detailed 
analysis of the Dyson- Schwinger equation for the electron inverse propagator which will lead to 
an examination of the eigenvalue problem in QED. The last section, Sec. 4 will be devoted to the 
conclusion that there is no non-trivial solution to the eigenvalue in finite QED followed by some 
important remarks. 



2 Dyson-Schwinger equation at the QED fixed point 

To investigate the Dyson-Schwinger equation of QED at the Gell-Mann-Low fixed point, we begin 
with the renormalized electron propagator whose general invariant form 

s-\p)=Mip') + np'), (2) 



follows from Lorentz invariance and parity invariance (PT invariance will suffice [|I^ ). In the finite 
QED theory, we set S(p^) = 0, for any p and chiral symmetry remains an exact sjniimetry. The 
proper renormalized vertex function in QED satisfies the Ward-Takahashi identity 

{p-p'n,{p,p') = s-\p)-s-\p'). (3) 

We would like to refer the reader to our earlier work [^] where we have reviewed the well-known 
consequences of the Gell-Mann-Low eigenvalue equation iplx) = 0, which may or may not have a 
non-trivial zero at xq = ao at the position of the bare fine structure constant of QED, ao = Z^^a. 
The important premise of the finite theory of QED is that the position of the zero, x = Xq can 
be determined by working with QED with zero physical mass [||. This is predicated upon the 
application of Weinberg's theorem [0, which ensures that terms vanishing asymptotically in each 
order of perturbation theory in the massive case do not sum to dominate over the asymptotic parts. 

It can be easily shown that at the Gell-Mann-Low fixed point with m = in finite QED, the 
full, exact, remormalized electron propagator has the simple scaling form WM 



s-\p)=M{p') = M^] (4) 




where 7(a) is the anomalous dimension of the electron in the massless theory given in the Landau 
gauge by 

'y = fi^\nZ2 = 0{a') + ---. (5) 

and fi is the subtraction point. This can be established as follows. If we begin with the Callan- 



Symanzik renormalization group equation [|15[, specialize to the Landau gauge and set m = 
(massless electron) and /3(a) = at the fixed point, then we have the equation satisfied by the 
two-point function which is essentially the same as the inverse electron propagator 

L|- + 27'|r(2)(p,a,0,/i,0) = 0. (6) 

The solution for the two-point function can be expressed as 

A(p^)^g)'. (7) 

This is customarily expressed in terms of the Euclidean momenta, as a function of (— p^), which is 
how we shall employ it later in this investigation. 

This can be confirmed by examining the trace anomaly in QED |l^. At a fixed point, (3{a) = 0, 



when we set the physical electron mass equal to zero, we have for the divergence of the scale current: 

d^D, = ^F.^F^'^ + [1 + 7e(a) ]m^V^ = (8) 

and hence scale invariance becomes exact. We further assume that scale invariance is not sponta- 
neously broken, Qd|0 >= 0, where 

Qd = I d^xDo{^,t), (9) 

{i.e., no dilatons are present in QED). Therefore Eq.(^) follows. 

We now consider the Dyson-Schwinger equation satisfied by the inverse of the full, exact, renor- 
malized electron propagator in m = QED: 

S-\p) = Z2i) - iZ2e\2ii)-^ J d^k YD^,S{k)T''{p, k). (10) 



It is well-known that in a theory of m = 0, spin-^ QED, the full, exact, renormalized photon 



propagator is exactly given by the free photon propagtor, as established by Eguchi |jT^ and thus 
we have the form for the photon propagator in the Landau gauge: 

where q^ = P^j, — k^. The solution to the Ward-Takahashi identity satisfied by the renormalized, 
proper vector vertex function in QED can be determined in the standard manner by the gauge 
technique [jl8| to yield 

r = r^ + r^, (12) 

where the longitudinal part of the vertex function admits the general, kinematical singularity-free 
solution [ H^, |20| given by 



mP,k) = l{A + A)Y + ^^,{l''^+hn 

- ^(^37)[2|^^'^^ + (P' + ^V], (13) 

and the transverse piece obeys the condition 

{p-p'),tUp,p') = o (14) 

and consequently undetermined by the Ward-Takahashi identity. Here we have employed the no- 
tation A = y4(fc^) etc. The transverse piece may be expressed by a kinematical-singularity free 
decomposition in terms of invariant functions but we find it is not necessary to do so at this point. 
However it must be stressed that we shall retain it throughout our calculation. We shall see that, 
undetermined and arbitrary as it is, the transverse piece of the vertex has a significant bearing on 
our conclusions. In massless finite QED, when the chiral symmetry is exact, the above solution for 
the vertex function reduces to 

r^(p, k)=^-{A + A)Y - \ (^^) [ 2|^7^^ + (/ + A:^)7^ ] + r^(p, k). (15) 

The Dyson-Schwinger equation, Eq. ([TO|) reduces to 



s-\p) = Mip') = U - iz,e\27r)-' J d'k [i^ - 7.) ^^ [ri(p, k) + r^ 



(p,fc)]. (16) 



We may now multiply the above equation by |^, divide by 4j9^ and evaluate the trace over the Dirac 
matrices. After a tedious computation, we arrive at the following equation: 



A{p' 



-Z2 = ^Z,e\2.rj d^k-^^ 



p'^k^A{k^) 
^(A + i) [2p'e - {p . k){p' + e)] + ^^^{P ■ k){p' - k'f 



(A + i)(^-^) 



A~A 



q2 q^i^k"^ — p"^) 



4{p-ky-{p-k){p' + k') +TUp,k\p-k)}, (17) 



where the last term, F;^ arises from the trace calculation of the term containing the transverse 
vertex piece. 

We observe that this is a non-linear integral equation satisfied by the function A(p^) and it is 
the exact consequence of the Dyson-Schwinger equation for the electron propagator in finite QED 
at the Gell-Mann-Low fixed point since we have not introduced any approximations and we have 
not discarded the transverse piece. 

To proceed further, we need to evaluate the various angular integrals appearing in the above 
equation. For this purpose, let us first transform to Euclidean momenta by the following transfor- 
mations: 

d^k -^ id^k = i f k^dkdVt; / -^ -/; A;^ -^ -k'^; p ■ k ^ -p ■ k. (18) 

Identifying the various angular integrals which occur as Ji, /2, ■ ■ ■ which are defined in Appendix A, 
we obtain 



A(V 



Z, = -Z,e^2.rjkHk-^ 



p'^k^A{-k^) 
(A + A) hp^Ph - {p^ + k^)h ]-{p'- k''){A - A)L 



- {A + A)h- 






4/3 - (p2 + e)h] + F^( V, -k\ ~p-k)\. (19) 



It is understood that all momenta are Euclidean, defined by Pe = V~P^? ^-B = \/—k^ in what 
follows. Henceforward we shall drop the subscript E for Euclidean momenta in order to avoid 



clutter. Performing the angular integrations and making use of the results in Appendix A |21|], we 
arrive at the following result 



A{p') - Z, 



kdk- 



1 



IGtt^ Jo p'^A{k'^ 

(A{p^)-A{k^) \2a\p^ 

[ (p2 - fc2) 



(20) 



.2\2 



Plil 



O^ 



Plil 



a 



+ rUp',e 



where r^{p'^,k'^) arises from the angular integral of the transvserse vertex part, and a = p^^/p-^ 
and p< = min{p,k}, p> = max{p,k} (see Appendix A ). This equation contains the essential 
ingredients of finite QED, with no approximations nor additional assumptions. Now in order to 
ascertain whether the Dyson-Schwinger equations of finite theoy of QED, as we have developed thus 
far, admit of a solution to the Gell-Mann-Low eigenvalue equation, we proceed as follows. We may 
check the self-consistency of the theory constructed in this manner, of massless QED at the fixed 
point by making the replacement 



A{p') 



p_ 



2 \ 7 



(21) 



in accordance with Eq.(^ where 7 = 7(a) is the QED anomalous dimension in the m = theory. 
After some algebra, we thus obtain 



p_ 



2\ T 



Z^e' 



+ 



IGtt^ jo 
16^ 



' kdkl[{pyky-l]^^'^' 



3k^ 



p^{p'^ — fc2) 



+ rUp',k') 



kdk\[{p^/k 



2w 



^Jp2p(p2_^2)+^T(P,fcJ 



(22) 



where T^ and T^ represent the contributions arising from the transverse piece vertex function. With 
a change in variables, s = p^, fc^ = sx, this can be rewritten in the form 



(s//i^ 



+ 



327r2 
327r2 



dx 
I dx 



x{l -3x){x-'' -1] 
x{x — 3){x^"' — 1] 



(1 



X] 



TUs.x) 



TUs,x) 



(23) 



where Ffn and Ff. are contributions arising from the transverse vertex piece. For general values of 



7, we can evaluate the integrals p2[ and we obtain the result 



7 p^ 

[s/fi^y - Z2 = ^{3F(l,3,4;l)-F(l,2,3;l) + F(l,2-7,3-7;l) 

- 3F(1, 3 - 7,4 - 7; 1) + 3F(1, -1, 0; 1) - F(l, -2, -1; 1) 

+ F(l,7-2,7-l;l)-3F(l,7-l,7;l) + r?.(s)} (24) 

in terms of the hypergeometric functions, where r^(s) arises from the integral of the contribution 
from the transverse vertex piece. 

If we evaluate Eq.(p^) at s = /i^, we obtain 

Z2 = l - ^{3F(l,3,4;l)-F(l,2,3;l) + F(l,2-7,3-7;l) 
- 3F(1, 3 - 7, 4 - 7; 1) + 3F(1, -1, 0; 1) - F(l, -2, -1; 1) 
+ F(l,7-2,7-l;l)-3F(l,7-l,7;l) + r^(/i2)}, (25) 

which can be rewritten as 



e 



2 



^2^ = 1 + 7T^{3F(l,3,4;l)-F(l,2,3;l) + F(l,2-7,3-7;l) 
327r^ 

- 3F(1, 3 - 7, 4 - 7; 1) + 3F(1, -1, 0; 1) - F(l, -2, -1; 1) 

+ F(l,7-2,7-l;l)-3F(l,7-l,7;l) + r5.(/x2)}, (26) 

This has been obtained in the Landau gauge to all orders in a. We are now ready to analyze the 
results contained in Eqs. (p^, p6D , a major consequence of the Dyson-Schwinger equations of finite 
QED with massless electron at a fixed point and draw conclusions. 



3 Conclusion and Summary 



Let us examine Eq.(p3D and determine what are the allowed values of 7, the anomalous dimension 
in massless QED. From Eqs. (|2^) and (^Sj) , we obtain 

{s/^?Y = l-^{3F(l,3,4;l)-F(l,2,3;l) + F(l,2-7,3-7;l) 

- 3F(1, 3 - 7, 4 - 7; 1) + 3F(1, -1, 0; 1) - F(l, -2, -1; 1) 
+ F(l, 7 - 2, 7 - 1; 1) - 3F(1, 7 - 1, 7; 1) + r^(;i2)| 

327r2 

- 3F(1, 3 - 7, 4 - 7; 1) + 3F(1, -1, 0; 1) - F(l, -2, -1; 1) 

+ F(l,7-2,7-l;l)-3F(l,7-l,7;l) + r^(s)}, (27) 



+ ;^ {3F(1, 3, 4; 1) - F(l, 2, 3; 1) + F(l, 2-7,3-7;!) 



which simplifies to 



{s|^^Y = l-^,{TW)-Tl{s)}. (28) 



We note that either 7 7^ or 7 = 0. Let us first consider the case 7 7^ 0. Since 7 is a finite 
function of a, the left hand side of Eq.(p8D is finite. On the other hand the right hand side is unity 
since, according to Eq . (^61) , Z2^ diverges, which means that Zg"^ is infinite for arbitrary non-zero 
values of 7. Eq.(P5|) cannot be satisfied unless 7 = 0. We are therefore forced to conclude / P5| / that 
7 = and thus 'j = is the only solution . Incidentally, one can show explicitly that the divergent 
terms of the longitudinal part in Eq.(PB|) do not cancel except when 7 = 0. (see Appendix B for 
details). With the choice 7 = 0, it then follows from Eq.(p6D that Zg"^ is finite and 



Z-' = Z-\f,') = I + ^^TUfi') (29) 



and hence we determine that r^(/i^) is finite. 

When 7 = 0, we observe from Eq.(^), that 

Z^e' 



327r2 



{r^(/i2)-r^(.)} = 0, (30) 



where Z2 is given by Eq. (|29|) and hence r^(s) is finite. The only way to obtain a non-trivial 
eigenvalue, e^ 7^ 0, is if Tj.{s) = Tj^lfi^). If the latter were valid, then this transverse vertex piece is 

8 



a constant, independent of s which amounts to an artificial and contrived constraint on the arbitrary 
and unknown transverse vertex part. We therefore conclude that the fine structure constant has 
only the trivial eigenvalue. 

Returning to our main result, we see that 7 = is based on the non-pertubative method 
stemming from the Dyson- Schwinger equations, and is thus valid to all orders in a in the Landau 
gauge. It therefore follows that for weak coupling {a << 1), a must vanish identically since 7 in the 
Landau gauge is of order a^. Next we may raise the question of what happens in strong coupling 
QED. Now we appeal to the gauge independence of the result. Since /?(«) is gauge independent 
(since Z3 is gauge independent) but 7(0;) is gauge dependent (since Z2 is gauge-dependent), the 
vanishing of /?(«) is valid in all gauges, but if 7 = in one gauge, it need not be so in other 
gauges. It then follows that the only gauge independent solution which is a simultaneous zero of 
both (3{a) = and 7(a) = is a = which is the trivial solution. We cannot invoke the Federbush- 
Johnson theorem [Q which only applies to the gauge independent sector |^ in QED. Since the 
gauge technique employed in this investigation is non-perturbative, we reiterate that our conclusion 



is valid for arbitrarily strong coupling p, |T0|, |TI 



We may recall that Adler et al |§] and Baker- Johnson [^ had argued from the triangle anomaly 
?7| that the Gell-Mann-Low function cannot have a zero at all. It is interesting that our 
present investigation establishes the absence of a non-trivial zero of the Gell-Mann-Low function 
by an entirely different approach [^ , in which the existence of the non- vanishing transverse vertex 
part plays a crucial role. 

We have carried out our investigation in the Landau gauge. The theory of QED is gauge 
invariant in content and we believe that our conclusions should prevail in any gauge. The form of 
the solution A{p'^) = (p^/fi^)'^ remains valid in all gauges in the minimal subtraction scheme ^9 



at /3(a) = 0. It would be interesting and worthwhile, however, to demonstrate the manifest gauge 
invariance of the theory. This complicated extension is under investigation and will be reported in 
a future communication. 



Appendix A: Angular integrals 



The computation of the right hand side of Eq.(|T9|) requires the evaluation of several angular 



integrals. A few but not all of these integrals are contained in |21|. We employ an expansion in 
terms of the Gegenbauer functions, 

1 °° 1 v 

where p> = max{p,k},p^ = min{p,k}. The evaluation of the angular integrals is facilitated by 
employing the standard properties of the Gegenbauer functions |3^, including in particular, the 
following: 

rl 



dxvT—x^C^{x)dx 



TT 



JnO, 



dxVl — x'^Cl^{x)C^{x)dx 






xC^{x) 



Cn+1 + ^n-li^) 



The integrals identified below can thus be evaluated and the results are as follows. 



h 
h 
h 
h 
h 



dQ 



\p — fcp 
p ■ kdQ 
\p — fcp 
{p ■ kfdQ 



2^- 
P> 



2 2 

71 a 



\p-k\^ 

dn 27r2 



2 2/1, 2\ 



Pi (1 - a^) 



\p-k\^ 

p-kdVL _ 27r2 (t2 

\p — k\^ Py 1 — cr^ 



(31) 



(32) 



The following notations have been used here: p> = max{p, /c},p< = min{p, k}, a = p</p>. It 
is understood that p, k stand for the Euclidean lengths pe = \/ —p^-, kE = V—k^ respectively. 



10 



Appendix B: Analysis of the hypergeometric functions 



Here we shall analyze the behavior of the hypergeometric functions appearing in Eg. (p6D . For 
this purpose we focus our attention on the coefficient of (e^/327r^) which is manifestly divergent on 
the right hand side of the equation. This may be expressed in terms of the gamma functions as 
follows: 

i mm _ r(3)r(6) r(3-7)r(6) _ T{A-j)m 
\ r(3)r(i) r(i)r(i)""r(2-7)r(i) r(3-7)r(i) 
rfe)r(e) 



+ 3 



^^-' + '^^^'^ + ^Jf^-^Ml - sjm^ + rW) 



r(i)r(-i + e) r(-2 + e)r(i) r(7-2)r(i) r(7-i)r(i) 

where lim e — > is understood. Each of the divergent terms can be dealt with by using the identity 



derived by Ryder |31| for the divergent gamma function 



where 



r(-n + e) 



%l){n 



-11 



n\ 



-+^{n+l) + 0{t) 
e 



-I 



1 

+ - 
n 



7e, 



is the logarithmic derivative of the gamma function [Q and 



7e 



lim 5Z Inn = 0.5772156649 



\m=l '"' 



is the Euler-Mascheroni constant. The conclusion stated in the beginning of Section 3 immediately 
follows, namely: Z^^ is infinite when 7 7^ 0; Z^^ is finite only if 7 = 0. Note that Vlp[jji^) is finite 
when 7 = 0. 
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